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Abstract 



We study the resummation of large logarithmic perturbative corrections to the partonic 
cross sections relevant for the process pp — > hX at high transverse momentum of 
the hadron h, when the initial protons are longitudinally polarized. We perform the 
resummation to next-to-leading logarithmic accuracy. We present numerical results 
for center-of-mass energies = 19.4 GeV, relevant for comparisons to data from 
the Fermilab E704 experiment, and \/S = 62 A GeV, where preliminary data from 
RHIC have recently become available. We find significant enhancements of the spin- 
dependent cross sections, but a decrease of the double-spin asymmetry for the process. 
This effect is less pronounced at the higher energy. 



1 Introduction 



The spin structure of the nucleon continues to be a particular focus of modern nuclear and par- 
ticle physics. As is well known, the total quark and anti-quark (summed over all flavors) spin 
contribution to the nucleon spin was found to be only about ~ 25%, implying that the gluon 
spin contribution and/or orbital angular momenta may play an important role. There is currently 
much experimental activity aiming at further unraveling the nucleon's spin structure. One em- 
phasis is on the determination of the spin- dependent gluon distribution, Agf, of the nucleon, which 
ultimately would give the gluon spin contribution to the nucleon spin. Deep-Inelastic scattering 
(DIS) has provided most of the presently available information on nucleon spin structure, but has 
left A^f essentially unconstrained [H El E]. Particularly good prospects for determining Ag{x, Q"^) 
over a wide range of momentum fractions x and scales Q are offered at the Relativistic Heavy-Ion 
Collider (RHIC) at BNL, which is the first polarized proton-proton collider. Spin asymmetries in 
high-energy pp scattering can be particularly sensitive to A(7, for processes where gluons in the 
initial state contribute already at the lowest order of perturbation theory ^ . One example is the 
single-inclusive production of large transverse-momentum [pt) hadrons, pp —>■ hX. Indeed, RHIC 
data taken at = 200 GeV on the double-spin asymmetry An for pp nX [5] and for the 



related process pp jetX [6] are now starting to put significant constraints on A^f, indicating 
that Ag is not too sizable in the accessed region of gluon momentum fractions. Similar conclusions 
are drawn from results obtained in lepton scattering [7j. 

RHIC is, however, not the first place where the spin asymmetry A'^j^ for pp — > nX was inves- 
tigated. The Fermilab E704 fixed-target experiment presented measurements of for 200 GeV 
protons impeding on a proton target [8J, resulting in a/S' = 19.4 GeV center-of-mass (cm.) en- 
ergy. An asymmetry consistent with zero was found for pions produced with transverse momenta 
1 < Pt ^ 4 GeV at central cm. system angles. An interesting question is whether this information 
already puts a constraint on Ag at the x values relevant here, 0.1 < x < 0.4. In [B] the experimen- 
tal data were also compared to theoretical leading-order (LO) calculations using various different 
Ag distributions. It was found that indeed there was some sensitivity of the data to Ag, with 
extremely large Ag (of size similar to the unpolarized gluon distribution in the accessed x region) 
seemingly ruled out. On the other hand, there are arguments against such a direct interpretation. 
For typical fixed-target kinematics as those in the E704 experiment, unpolarized single- inclusive 
hadron cross section data are generally not at all described even by next-to-leading order (NLO) 
(let alone, LO) theoretical calculations [9], with theory falling way short. One may therefore won- 
der if it is then sensible to confront LO calculations for with the data. In order to address 
this issue, Ref. [10] considered the effects of possible Gaussian-distributed "intrinsic" transverse 
momenta (fc-r) of the partons on A'J^^ in this kinematic regime. It was found that intrinsic 
tends to decrease the spin asymmetry significantly, so that even relatively large Ag appeared to 
be compatible with the E704 data. At the same time, intrinsic /c^ improves the comparison with 
the unpolarized cross section data. However, from a theoretical point of view, implementation 
of intrinsic kx into a single-inclusive cross section is not really a satisfactory approach because 
one only has collinear factorization in this case. At best, intrinsic kx effects may be regarded as 
providing an effective model for possible power-suppressed contributions to the cross section. Im- 
plementation of intrinsic k^ also obscures the role of perturbative higher-order contributions to the 
cross section. Nonetheless, the results of [10] indicate that there can be substantial contributions 
to Ali^ in the fixed-target regime that go beyond low orders of perturbation theory. 
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Progress on the theoretical description of the unpolarized single-inclusive hadron cross section 
in the fixed-target energy regime was made in Ref. |Tl]. For typical fixed-target kinematics, 
the value of xt = 2pt/VS is relatively large, xt ^ 0.1. It turns out that the partonic hard- 
scattering cross sections relevant for pp —>■ hX are then largely probed in the "threshold" -regime, 
where the initial partons have just enough energy to produce the high-transverse momentum 
parton that subsequently fragments into the hadron, and its recoiling counterpart. Relatively little 
phase space is then available for additional radiation of partons. In particular, gluon radiation is 
inhibited and mostly constrained to the emission of soft and/or coUinear gluons. The cancellation 
of infrared singularities between real and virtual diagrams then leaves behind large double- and 
single-logarithmic corrections to the partonic cross sections. These logarithms appear for the first 
time in the NLO expressions for the partonic cross sections, where they arise as terms of the form 
as ln^(l — x^) in the rapidity-integrated cross section, where xt = Pt/V^ with px the transverse 
momentum of the produced parton and s the cm. energy of the initial partons. At yet higher 
(/cth) order of perturbation theory, the double-logarithms are of the form a|ln^'^(l — x^). When 
the threshold regime dominates, it is essential to take into account the large logarithms to all 
orders in the strong coupling as, a technique known as "threshold resummation" |T2]. Based 
on earlier work [131 IE] oii the resummation for 2 — >■ 2 QCD hard-scattering, we examined the 
effects of threshold resummation on the single-inclusive hadron cross section in [11] and found 
very significant enhancements of the theoretical prediction in the fixed-target regime, which in 
fact lead to a relatively good agreement between resummed theory and the data. This also sheds 
light on the size of additional power-suppressed contributions to the cross sections (among them, 
perhaps, effects related to intrinsic kx] see also Ref. [15]), which do not seem to play a dominant 
role. We concluded that threshold resummation is an essential part of the theoretical description 
in the typical fixed-target kinematic regime. Its effects at higher energies (such as at RHIC) are 
much smaller, even though it has to be said that one is here typically much further away from the 
threshold regime so that the applicability of threshold resummation is not entirely clear then. 

In the light of the results of Ref. [H] and the E704 data, it appears desirable to apply threshold 
resummation also to the spin asymmetry A^^, which is the goal of this paper. In this way, one 
may hope to put the theoretical description of A'^^ for single-inclusive hadron production in the 
fixed-target regime on firmer ground. One may then also revisit the question as to whether the 
E704 data already allow to put a constraint on Ag. 

We also note that recently preliminary data for the cross section and double-spin asymmetry 
taken at RHIC's lower energy = 62.4 GeV have been reported [HI HTj. Even though the 
approximations needed for threshold resummation to be useful work slightly worse for the kine- 
matics relevant here, it is of great interest to confront the resummation with the data. This will 
also be done in this paper. 

The remainder of this paper is organized as follows: Section [2] summarizes the theoretical 
perturbative-QCD framework for the process under study. In Sec. [3l we present the resummed 
spin-dependent cross section to next-to-leading logarithmic (NLL) accuracy. Section HI gives phe- 
nomenological results for the effects of threshold resummation on the spin-dependent high-p^ pion 
cross section at \/S = 19.4 GeV and at a/S* = 62.4 GeV, and on the corresponding double-spin 
asymmetries Aj^j^. Finally we draw our conclusions in Sec. O Two Appendices contains the 
relevant ingredients for the resummation in the spin- dependent case. 
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2 Cross section and spin asymmetry in perturbation 
theory 



We are considering the process 

p{p,, Ai) + Pip,, A2) Hps) + X , (2.1) 

where the Aj denote the hehcities of the initial protons, and the pi [i = 1,2,3) are the four- 
momenta of the "observed" hadrons. One defines the spin-averaged and spin- dependent cross 
sections as 

da = ^ [da{Ai = +, A2 = +) + da{Ai = +, A2 = -)] , 
dAa = ^ [da{Ai = +, A2 = +) - rfa(Ai = +, A2 = -)] , (2.2) 
respectively, and their double-spin asymmetry as 

dAa 

All = -^. (2.3) 
aa 

Hadron h is assumed to be produced at large transverse momentum pt- For such a large- 
momentum-transfer reaction, the factorization theorem [18] states that cross section may be fac- 
torized in terms of collinear convolutions of parton distribution functions for the initial protons, a 
fragmentation function for the final-state hadron, and short-distance parts that describe the hard 
interactions of the partons and are amenable to QCD perturbation theory. The long-distance par- 
ton distributions and fragmentation functions are universal, i.e., they are the same in any inelastic 
reaction. Long- and short-distance contributions are separated by a factorization scale. 



As discussed in Ref. [TT], a major simplification of the resummation formalism occurs when 
the cross section is integrated over all pseudo-rapidities t] of the produced pion. This will also be 
done in this paper. The factorized spin-dependent cross section for pp —>■ hX can then be written 
as 



p^dAa^xx) y-v 

a,b,c 



dpi 



dxi Afa {xi, iJ,"^) / dx2 Afb {x2, 1^^) / dzz^ Dh/c{z,^i^) 
Jo Jo 



Here the A/a h are the spin- dependent parton distributions in the proton, 

A/,(x,^2) = f+{x,f,')-f-{x,^,') (2.5) 

with /+ (/~) denoting the distribution of parton type a with positive (negative) helicity in a 
proton of positive helicity. Dh/c is the fragmentation function for parton c fragmenting into the 
observed high-p^ hadron. The sum in Eq. (12. 4p runs over all partonic channels, with the associated 
spin-dependent partonic cross sections dAaab^cX- The latter are defined analogously to Eq. (12. 2p . 
with helicities now corresponding to parton ones. They are perturbative and have an expansion 
of the form 

dAa^b^^x = dAa^l^^ + -rfAai;LcX + • • • (2-6) 

71 
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with as the strong couphng. The scale fi in Eq. (12.41) is the factorization scale. We could distinguish 
in principle between factorization scales for the initial state (parton distributions) and the final 
state (fragmentation function). For simplicity, we will not do this in this paper. There is also 
a renormalization scale, at which the strong coupling constant is evaluated. We will collectively 
denote all scales by /i. The dependence on /i is implicit in the partonic cross sections in Eq. (12. 4p . 
Finally, fj is the pion's pseudorapidity at parton level, related to the one at hadron level by 



fj = rj — ^ ln(a;i/a;2). Its limits are given by r)+ = —fj^ = In (1 + a/1 — x\)/xt 
xt = 2pt/VS, and its partonic counterpart is xt = 2p^/\^ = xt j z^x^x^- 



where, as before. 



We note that the corresponding expression for the factorized spin-averaged cross section is 
obtained from Eq. (12. 4p by dropping all A's, meaning that the spin-dependent parton distributions 
are replaced by their usual unpolarized counterparts, and the partonic scattering cross sections 
by the spin-averaged ones. 



3 Resummed cross section 



As mentioned above, we will follow |TT] to perform the threshold resummation only for the case 
of the fully rapidity-integrated cross section. The resummation of the soft-gluon contributions is 
achieved by taking a Mellin transform of the cross section in the scaling variable x\.: 

AaiN)^ Cdxlixl)''-' Pl^P^. (3.7) 
Jo dpT 

For the rapidity-integrated cross section, the convolutions in Eq. (12.41) between parton distri- 
butions, fragmentation functions, and subprocess cross sections then become ordinary prod- 
ucts [II1[I9]: 

Aa{N) = ^ AUN + 1, fi^) AMN + 1, fi^) Dh/,{2N + 3, fi^) Aaab^cxiN) , (3.8) 

a,b,c 

where 

Aaab^cx{N) = dxrr [xt) J dr] — . (3.9) 

In Mellin-moment space, the threshold logarithms become logarithms in the moment variable A^. 
The leading logarithms are of the form \n^^ N; subleading ones are down by one or more powers 
of IniV. Threshold resummation results in exponentiation of the soft-gluon corrections in moment 
space p!2| [13] . The leading logarithms are contained in radiative factors for the initial and final 
partons. Because of color interferences and correlations in large-angle soft-gluon emission at NLL, 
for QCD hard-scattering the resummed cross section becomes a sum of exponentials, rather than 
a single one [HI , unlike the much simpler cases of the Drell-Yan or Higgs cross sections ^12j . 

Combining results of [121 [131 [20], we can cast the resummed spin-dependent partonic cross 
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section for each subprocess into a relatively simple form [11^ 



TT^ T)^ TT^ T'^ 



ab^cd N 



(3.10) 



where Aa^^^^l{N) denotes the LO term in the perturbative expansion of Eq. (13.91) for each process. 
We list the moment-space expressions for all the spin-dependent Born cross sections in Appendix 
A. Each of the functions Jff,V%,vff''''^'"^ in Eq. (^M) is an exponential. represents the 
effects of soft-gluon radiation collinear to initial parton a and is given, in the MS scheme, by 

r.(l-^)2Q2 ^ 

~i2 



1 ^N-l 



1 







z 



^Aa{asiq'^))dz 



(3.11) 



and similarly for A^. Here, = 2p^. We will specify the function Aa below. Collinear soft-gluon 
radiation to parton c yields the same function [20]. The function Jf^ embodies collinear, soft or 
hard, emission by the non-observed recoiling parton d and reads: 



In 4 



1 ^N-l 



1 



{1-^)Q2 



dq 



1 



{i-^)2Q2 q 



^AMs{q')) + ^BMs{{l-z)Q')) 



dz . 



(3.12) 



Large-angle soft-gluon emission is accounted for by the factors which depend on the 

color configuration / of the participating partons. A sum over the latter occurs in Eq. (13.101) . with 
AGab^cd representing a weight for each color configuration, such that ^/ AG^^,^^^ = 1. Each of 
the U)^ IS given as 



In -T-)(int)afe^cc( 
ni 1^1 



1 



-Diab^cdiasi{l-zYQ^))dz . 



(3.13) 



Finally, the coefficients ACab~*cd contain A^— independent hard contributions arising from one-loop 
virtual corrections. 



Most of the ingredients to Eqs. (I3.11l) - (l3.13p are well-known from the literature because they co- 
incide with the results obtained for spin-averaged scattering. This is the case for the functions Aa, 
Ba, and Djab^cd, because these are associated with soft gluon emission, which is spin-independent. 
The only differences between the spin-dependent and the spin-averaged cases reside in the coef- 
ficients AG^^_^^^, ACab^cd and of course in the Born cross sections. These terms are all related 
to hard radiation, which depends on the polarization state and therefore in general differs for the 
polarized and unpolarized cases. 



We first briefiy recall the known functions and then turn to the new parts. Each of the functions 
J-' = Aa, Ba, Diab^cd is a perturbative series in as. 



with ED: 



J^(as) 



TT V TT / 



as 



= c 



Ca 



Ca 



67 
18 



TT 

"6 



Bi'^ 



la 



(3.14) 



(3.15) 



^Note that the symbols 'D)^,'Dj 



(int)afc — '^cd 



N 



in the equation below are usually referred to as A^vjA^ 



(int)a6 — 



N 



the literature [TT]. We have changed this notation in order to avoid confusion with the label "A" indicating 
spin-dependent cross sections and parton distributions in this paper. 
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where Nf is the number of flavors, and 

Cg = CA = N, = 3, C, = Cf = {Nl - l)/2iV, = 4/3 

7, = -?>Cf/2 = -2 , 73 = -27r6o , h = {HCa - 2Nf) . (3.16) 

The coefficients D^]^lh^cd hsted in Ref. [TT] . 

In order to determine the effects of the color interferences for large-angle soft-gluon emission 
in the 2^2 processes ab cd, we can follow the procedures presented in [THl [S] . The soft- 
anomalous dimensions and soft factors determined in [13] are again identical in the spin-averaged 
and the spin-dependent cases. The differences arise solely in the color-connected Born cross 
sections. We therefore only need to derive the latter for polarized scattering, using the same color 
basis as that chosen in [13] . The results in [13] have actually been given for arbitrary rapidity; for 
the case of the rapidity-integrated cross section we consider here it is sufficient to set 17 = (see [TT] 
for more detail). Nonetheless, for future convenience, we present in this work the spin- dependent 
color-connected Born cross sections also at arbitrary rapidity. In this way, they may be directly 
used in future investigations of the resummed cross section at fixed rapidity. The results are 
collected in Appendix B. For the case of the rapidity-integrated cross section, the color-connected 
Born cross sections at 17 = 0, when normalized to the full Born cross section for each par tonic 
channel, give the color weights AG^^^^^ [TT] that we need to NLL, which are listed in Appendix A. 

The perturbative expansion of the coefficients ACab^cd reads: 

Aa,_,. = 1 + ^ACiJLe. + 0{al) . (3.17) 

TT 

In order to determine the coefficients AC^),^^^, we take advantage of the full analytic NLO calcula- 
tion of Ref. [22] . For each partonic channel one expands the resummed cross section in Eq. (13.101) 
to ffist order in as. Near threshold, one can straightforwardly take Mellin moments of the full 
NLO expressions of [22]. By comparison of the two results one first verifies that all logarithmic 
terms in the full NLO results are correctly reproduced by the resummation formalism. The re- 
maining A^-independent terms in the NLO cross section give the coefficients AC^)^^^^. These turn 
out to have rather lengthy expressions, and we only give them in numerical form in Appendix A. 

This completes the collection of the ingredients for the resummed partonic cross sections. In 
the exponents, the large logarithms in occur only as single logarithms, of the form In'^'^^(A^) 
for the leading terms. Next-to- leading logarithms are of the form agln'^(A^). Knowledge of the 
coefficients given above allows to resum the full LL and NLL full towers in the exponent. It is 
useful to expand the resummed exponents to definite logarithmic order: 

lnP^(as(/i'),QV/i') = lnNh^^\X) + h'^^\X,Qyfi^)+0{as{aslnN)') , (3.18) 
In JUas{fi'),QVfi') = InN ff\X)+ff\X,Qyf,') + 0{as{aslnNf) , (3.19) 

where A = boasifJ^'^) In N . The functions /i*^*-* and /'•*•' have for example been given in [TT]. h^^^ 
and Z^^-* contain all LL terms in the perturbative series, while /i^^-* and f^'^^ are only of NLL 
accuracy. For a complete NLL resummation one also needs the coefficients In "D^ whose 
NLL expansion reads: 

r)(i) 

lnPK^"(as(/.?,),QV/^«) = ln(l - 2X) + O (as(as In Ar)'=) . (3.20) 
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In order to obtain a resummed cross section in space, one needs an inverse Mellin trans- 
form. Here one has to deal with the singularity in the perturbative strong coupling constant in 
Eqs. f l3.lip - fl3.13p . which manifests itself also in the singularities of the functions /i*^^'^) and /''^'^■' 
above at A = 1/2 and A = 1. We use the Minimal Prescription developed in Ref. |23], which 
relies on use of the NLL expanded forms Eqs. f l3.18p - p.20p . and on choosing a Mellin contour in 
complex- space that lies to the left of the poles at A = 1/2 and A = 1 in the Mellin integrand: 



dpT 



dN 

27Ti 



(3.21) 



where 6o«s(/^^) IhCa/p < 1/2, but all other poles in the integrand are as usual to the left of the 
contour. The result defined by the Minimal Prescription has the property that its perturbative 
expansion is an asymptotic series that has no factorial divergence and therefore no "built-in" 
power-like ambiguities. 



Finally, in order to make full use of the available fixed-order cross section [221 [Mj . which in 
our case is NLO (0(ag)), we match the resummed cross section to the NLO one. We expand the 
resummed cross section to 0{al), subtract the expanded result from the resummed one, and add 
the full NLO cross section: 



dpT 



Cmp+ioo 



2m 



(4) 



X 



(res) 



*cd 



(AT)- A 



(res) 
ab^cd 



{N) 



+ 



dpT 



(3.22) 



where A(3"^^^^^^(A^) is the polarized resummed cross section for the partonic channel ah — >■ cd as 
given in Eq. fl3.10p . In this way, NLO is taken into account in full, and the soft-gluon contributions 
beyond NLO are resummed to NLL. Any double-counting of perturbative orders is avoided. 



4 Phenomenological Results 

We are now in the position to present numerical results for the threshold-resummed spin-dependent 
cross section and spin asymmetry in single-inclusive pion production in hadronic collisions. We 
will focus here on pp — > tt^X in fixed-target scattering at \fS = 19.4 GeV and at = 62.4 GeV 
at the RHIC collider. In both these cases, experimental data exist [8| [T7]. 

For our calculations, we need to choose sets of parton distribution and pion fragmentation func- 
tions. To study the sensitivity of the measured spin asymmetries to the spin-dependent parton 
distribution functions, in particular the gluon density Ag, we use the "Gliick-Reya-Stratmann- 
Vogelsang (GRSV)" pLj and the "de Florian-Sassot (DS)" [2] densities. These both offer various 
sets of distributions, distinguished mostly by Ag distributions of different sizes. The set labeled 
"GRSV" is the regular ("standard") GRSV set. We will also use the GRSV "max G" set, which 
has a much larger Ag, given by assuming Ag{x) = g{x) at the initial scale for the parton evo- 
lution. Likewise, the DS "i-f-" was constrained to have a much lower Ag than the "iii-|-" one. 
For the spin-averaged cross section, we employ the MRST2002 set throughout. The pion 
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10^ 
10-6 

2 3 4 5 6 

p^(Gev) 

Figure 1: NLO and NLL resummed cross sections for polarized pp ^ tc Ox 75 = 19.4 GeV,for 
various sets of spin- dependent parton distributions of |21 IMI- We also show the 0{al) expansions 
of the resummed cross sections, and the analogous results in the unpolarized case. For better 
visibility, we have applied numerical factors to some results, as indicated in the figure. In the 
upper right inset, we present the ratios between the NLL resummed cross sections and the NLO 
ones. 

fragmentation functions are taken from the most recent analysis of e^e~ and pp data, "de Florian- 
Sassot-Stratmann (fDSS)" Ref. [26]. It is worth noticing that, according to Eq. (13.221) . one would 
like to have the parton densities and fragmentation functions available in Mellin-moment space. 
Technically, since most of the distributions are only available in x space, we first perform a fit with 
a simple functional form to the distributions, of which we are able to take moments analytically. 
This has to be done separately for each parton type and at each factorization scale. 

In Fig. [1] we present our results for the spin-dependent cross section at = 19.4 GeV, 
integrated over all pion pseudo-rapidities 77, for the various sets of the polarized parton distributions 
of PQ [2]. We have chosen all scales as /i = pt- We show the full NLO cross section based on the 
calculations in P2], as well as the NLL resummed predictions. We also display the expansions of 
the resummed cross section to 0{a\), which is the first order beyond LO. As can be observed, 
these faithfully reproduce the full NLO result, implying that the threshold logarithms addressed 
by resummation indeed dominate the cross section in this kinematic regime, so that resummation 
is expected to be useful. Towards lower ~ 3 GeV, the expansions slightly overestimate the 
NLO cross section, which is expected since one is further away from the threshold regime here, so 
that the soft-gluon approximation tends to become less reliable. For the sake of completeness the 
unpolarized case is also represented in Fig. [H 

The inset in Fig. [1] shows the resummed "i^'-factors" for the cross sections, defined as the ratios 
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1 2 3 4 5 6 7 

p^(Gev) 



Figure 2: Same as Fig. [I] hut for \fS = 62.4 GeV. 
of the resummed cross sections to the NLO cross sections (polarized or unpolarized) : 

As can be seen, 

^(res) 

is very large, meaning that resummation results in a dramatic enhancement 
over NLO. For the unpolarized case, this finding is in line with our previous results in |TT]. It is 
interesting to see that K^^'^^^ is large also for all sets of spin-dependent parton distributions. It is 
evident, however, that the enhancement is somewhat smaller than in the unpolarized case. This 
immediately implies that the spin asymmetry will generally be reduced when going from NLO 
to the NLL resummed case. One also notices that the resummation effects vary slightly for the 
various polarized parton densities. This may be understood from the fact that A^f is of different 
size in the various sets. Typically, resummation effects are more important for partonic channels 
with more external gluons [TT], so the size of A^f matters. 

Fig. [2] shows similar results for pp vr^X at = 62.4 GeV. As expected from the fact that 
one is further away from threshold here, the soft-gluon approximation becomes somewhat less 
accurate in this case, in particular at the lower px- One can observe that the resummation effects 
are generally much smaller at \^ = 62.4 GeV than at = 19.4 GeV. 

As mentioned earlier, we have determined the resummed formulas for the fully rapidity- 
integrated cross section, whereas in experiments typically only a certain limited range in rj is 
covered. In order to be able to compare to data, we therefore approximate the cross section 
(polarized or unpolarized) in the experimentally accessible rapidity region by 

p3^^{match) _ 3^^^(NLO) _ 

; (r] mexp. range) = il ^ ' ; (r/ mexp. range) , (4-24) 

dpT dpT 
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Figure 3: Results for the double-spin asymmetry at NLO and for the NLL resummed case for 
various sets of polarized parton distributions, at \/S = 19.4 GeV. We also show the experimental 
data of ig/. 

where K^"^*^^^ is as defined in Eq. 04.231] in terms of cross sections integrated over tlie full region 
of rapidity. In other words, we rescale the matched resummed result by the ratio of NLO cross 
sections integrated over the experimentally relevant rapidity region or over all rj, respectively. [27] 

Figure [3] shows our results for the spin asymmetry at = 19.4 GeV, for the NLO and 
NLL resummed cases, defined as in Eq. (12. 3p . averaged over the pion's Feynman-Xi? = XTsinh(?7), 
\xf\ < 0.1. Again the scales have been chosen to be /i = pt- We also show the data by the 
Fermilab E704 experiment [B]. As expected from Fig. [H generally decreases significantly 
from NLO to NLL. After NLL resummation, even a set with a very large Ag, such as the GRSV 
"maximal" scenario (which is now already ruled out by other measurements [5l[6l[7]) shows rough 
agreement with the data, given the rather large experimental uncertainties. It is interesting to 
note that similar results were found in |in] on the basis of LO studies invoking "intrinsic-fcr" 
effects. 

We now return to the case of pp scattering at RHIC at = 62.4 GeV. Recently, first 
preliminary data for the spin-averaged high-pr pion cross section as well as for the spin asymmetry 
^LL were reported by the Phenix collaboration [16j. The data cover the pseudo-rapidity region 
Ir^l < 0.35. Fig. H] compares our NLO and NLL resummed results for the spin-averaged cross 
section to the Phenix data. We use the scales /i = Cpt with ( = 1/2, 1, 2. It is interesting to see 
that the data lie at the upper end of the rather wide NLO scale band, whereas the resummed 
predictions have a smaller scale dependence and tend to describe the data rather well with scale 
II = Pt- We remind the reader that at y/S = 200 GeV the RHIC data are very well described by 
NLO with scale n = pt |5], while in the fixed-target regime resummation effects were found to 
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Figure 4: Invariant cross section for pp n at y/S = 62 A GeV at NLO and for the NLL 
resummed case. We also show the preliminary Phenix data /i 



be very significant (see [TT] and Fig. [T] above). We interpret all these features as indicating that 
threshold logarithms start to become relevant at = 62.4 GeV, which is "half way" between 
the typical fixed-target regime and RHIC's 200 GeV. Encouraged by the results in Fig. HJ we 
show in Fig. [5] our NLO and NLL results for the spin asymmetry Aj^. = 62.4 GeV, along 

with the Phenix data [17J. As can be seen in Fig. [5l sets with a large gluon polarization, like 
"GRSV max G", show a clear disagreement with the preliminary data. We observe that there 
is again a decrease of A\j^ when going from NLO to NLL, but that the resummation effects are 
somewhat smaller than what we found at \/~S = 19.4 GeV. We have explicity checked that very 
similar results for the asymmetries are obtained when implementing other sets of fragmentation 
functions, like those from [281 [29]. 



5 Conclusions 



We have studied in this paper the NLL resummation of threshold logarithms in the partonic 
cross sections relevant for the process pp —>■ hX at high transverse momentum of the hadron 
h, when the initial protons are longitudinally polarized. We have found that, like for the spin- 
averaged case pni, the resummation effects are large for the spin-dependent cross section in the 
typical fixed-target regime at ~ 20 GeV. The spin asymmetry Al^ is significantly reduced by 
resummation. A phenomenological consequence is that the Fermilab E704 data [8] are compatible 
with essentially all currently sets of spin- dependent parton distribution functions, among them 
sets with a rather large gluon polarization. 

We have also applied the resummation to the case \^ ~ 62.4 GeV, at which recently pre- 
liminary data from RHIC have become available for both the spin-averaged cross section and the 
double-spin asymmetry P^[T7] . For this case sets with a very large gluon polarization in the range 
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Figure 5: Results for the double-spin asymmetry at NLO and for the NLL resummed case for 
various sets of polarized parton distributions, at ^/S = 62.4 GeV. We also show the preliminary 
experimental data of fT7\l . 

0.05 ^ X < 0.2 are ruled out and distributions with a moderate to small gluon polarization are 
favored by the data. We find that resummation tends to lead to an improved description of the 
cross section. Its effect on the spin asymmetry is rather modest. We remind the reader that the 
threshold resummation is somewhat less reliable in this kinematic regime since one is further away 
from threshold than for the fixed-target case, so that subleading perturbative corrections may be 
more relevant. It will be desirable in the future to further improve the resummed calculation at 
\/S = 62.4 GeV by including terms that are subleading near threshold. More reliable conclusions 
regarding the effects of threshold resummation on A'J^^^ at this energy should then become possible. 
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Appendix 



A Results for the various subprocesses 

In this appendix we compile the moment-space expressions for the spin-dependent Born cross 
sections for the various partonic subprocesses, and the polarized process-dependent coefficients 
^^ab^cd^ ^C^iab-^cd Contributing to Eq. fl3.10p . As we mentioned in the main text, the coefficients 
ah-^cd are the same as in the unpolarized case; they have been given in Ref. [llj . Since the 
^ClfeLcd have rather lengthy expressions, we only give their numerical values for Nj = 5 and 
the factorization and renormalization scales set to jj, = Q. In all expressions below, Ca = 3 and 
Cf = {C\ - 1)/2Ca = 4/3. 



qq' qq': 



qq' qq': 



qq — i> q'q': 



qq qq: 



qq qq: 



^Giqqi^qqi = Giqq/^qq/, AG'2 gg'^gg' = G2qqi^qq' , 

^C^Z^qq, = 17.9311 {Nj = 5) . (A.l^ 



A-:?ri,(iV) = (SAT^ + SAT) B (at, , 

AGiqqi^q-/ = Giqqi^q-/, AG2qq'^qq' = G2qqi^qq' , 

^C^^qq' = 20.7021 {Nj = 5) . (A.2) 



^^(Born) = -a(^°''",\(Ar), 
qq~*q'q' ^ ' qq—*q'q' ^ ' ' 

AGJ^). = gS'1 ,-„ Ad'l ,-, =Cf^- (A.3) 

lqq—>q q 1 qq^q'q' ^ 1 99~><? <? ^IQ^Q Q 



A^K(iV) = {Ga{SN' + 5N) - 2N{3 + 2N)) B (^N, ^- 

AG,qq^qq = 7/5, AG 2 qq^qq = '2 / 5 , Acj^^^^^^ = 14.5364 (A^^ = 5) . (A.4) 

^^SZqliN) = «^^iV {Ga{2 + N){11 + 5N) - (A^ + 3)(5 + 2A^)) B (^N, , 
AGiqq^qq = "3/13 , AG2qq^qq = 16/13 , Acf'^^-^- = 24.2465 (A^^ = 5) . (A.5) 
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qq gg: 



A "(Born)/ _ _-(Born)/^^ .^{1) _ ^(1) 

^'^qg-^ggV^ ) - ^qq^ggV^ )^ ^'^Iqq^gg - '-'Iqq^ggy 

AGigq^gg = G I g g , AG2gq^gg = G 2 qq^ g g . (A. 6) 



qg qg: 



AGiqg^qg = Giqg^qg, AG2qg^qg = G2qg^qg, AG^qg^qg = G^qg^qg , 

AGil^qg = 14-2048 {Nf = 5) . (A.7) 



qg gq: 



qg^gq^ > 
AGi 



a 



+ 2^-4) (3A^' + ^N)B f iV, \ 



qg^gg 



Gil 



AGI 



(1) 



qg^gq 



-qg^gg^ 
21.2354 (A^/ 



AG2qg^gq 

5). 



"iqg^gq^ 



AG 



3 gg-^gg 



3 gg^gg ? 



(A. 



^ gg-- 



^^Igg^gg — ^Igg^ggj AG2 gg~tgg = G2 gg^gg, AG^gg^gg = G^gg^gg , 

^^igg-^gg = 20-3233 (A^^ = 5) . (A.9) 



gg QQ- 



^^gg^qqy^^ ) - ^gg^qqV^)i ^'-'Igg^qq- '-'Igg^qqi 

AGlgg^qq = Glgg^qq, AG2gg^qq = G2gg^qq. (A. 10) 



In the above expressions, B{a,b) is the Euler Beta-function. 



B Spin-dependent color-connected Born cross sections 



In this appendix we compile the color-connected Born cross sections that we need for our study. 
Our choices for the color bases follow precisely [13]. For a given color basis, the color-connected 
Born cross sections will appear as a matrix that we shall refer to as "hard matrix". We will 
not repeat the expressions for the soft matrices S, and the anomalous dimension matrices F 
in these bases, which may all be found in [13]. Like [l3j, we present our results for arbitrary 
partonic rapidity, even though for our actual study of the rapidity-integrated cross section we 
only need the case 17 = 0. For each partonic reaction ab —>■ cd, we define the Mandelstam variables 
■s = (Pa + PbY, t = {Pa — Pc)"^, u = {pa — PdY ■ t,u 3X6 fuuctious of fj. In all expressious below, 
iV, = 3 and Gf = {N^ - l)/2N^ = 4/3. 
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We begin with the quark-antiquark annihilation processes. Depending on the quark flavor, 
there are three different quark-antiquark subprocesses to consider, q^qj — > qjqj , qjqj — > q^qk and 
QjQk QjQk- Each of these have their own hard matrix elements: 



--"12 



st 



(B.l) 



where H°'j^'^'^ refers in each case to the corresponding hard matrix for the unpolarized case given 



in m. 



Here one has simply: 



(B.2) 



QjQk QjQk 

The polarized hard matrix at lowest order can be expressed as. 



Qk 



a. 





2{s^ -u^)/{N^e) 



(B.3) 



There are only two different quark-quark processes to consider, depending on the quark flavors: 
QjQk QjQk 

Here the result is the same as for the process qjqt — ^ IjQk, but with the color basis inverted 
(equivalent to an interchange of the entries in the first and second rows and columns, Ci C2). 



qjqj q^qj 
Here, 



^-"12 



a. 



a 



a 



iVc-- 

tu 



'iV2 

c 

1^ 



+ 



1 is^-e) 



NrtU 



V? 



(B.4) 
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In the spin-dependent case one has: 



2{s^ -u')/{Nie) 




(B.5) 



99 

Here, 



qq-^gg 



qq^gg 



(B.6) 



99 m 

Again, the polarized coefficients are the negatives of the unpolarized ones: 

^j^gg^qq _ _]^gg^qq 



(B.7) 



(19 (19 
Here, 



c2 — ^,2 

ATj-qg^qg _ " ugg^qg 



Finally, we consider gluon-gluon scattering. For simplicity, we set = 3 explicitly here. 

99 99 

e4 _ +4 _ 4 

Afrgg^gg _ rrgg-^qg 
s^ + t^ + u^ 



(B.8) 



(B.9) 
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